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1. Introduction 

Chiral phase transition at finite temperature T c is one of most prominent features of QCD. In 
1980’s, Psarski and Wilczek investigated the order of the phase transition applying the e expansion 
to linear sigma model (LSM) |[TJ], and pointed out that the nature of the chiral phase transition of 
massless two-flavor QCD depends on the existence of the flavor-singlet axial symmetry, £4(1), at 
the critical temperature. 

Two extreme cases have been well studied so far. If the £4(1) breaking is infinitely large at 
T c , the system would be described by the standard 0(4) LSM and thus undergo second order phase 
transition. On the other hand, if the breaking vanishes at 74 [j2[] , the symmetry of the system is 
enhanced to U (2) x U (2). In this case, the nature of the transition is still under debate. It is argued 
that if it is of second the critical exponents in the £/(2) x U (2) universality class will emerge ii- 

In this work, we investigate chiral phase transition of massless two-flavor QCD with a finite 
£4(1) breaking. Instead of dealing with two-flavor QCD directly, £/(2) x U (2) LSM with a £4(1) 
breaking, which we call the £4(1) broken LSM below, is analyzed as the corresponding effective 
theory. The nature is investigated through the renormalization group (RG) flow calculated in the 
framework of the e expansion. Since the full results on this study are available in Ref. [ 4 1 , we focus 
only on the highlights. In the study of phase transitions, it is common to see effective potential. We 
present the preliminary result of the effective potential. 

2. RG flow of the £4(1) broken LSM 

The LSM Lagrangian we study consists of the part preserving £/(2) x U (2) symmetry and the 
one breaking £4(1), 

= -^U(2)xU(2) T -^breaking, (2-1) 

■%’u( 2 )xu( 2 ) = ^tr[^d><9 p <L T ] + ^m 0 2 tr[<Lfl) t ] + yg, (tr[<M> f ]) 2 + yg 2 tr [(<$><L» V ) 2 ] , (2.2) 

Ca 71^ 71^ 

•^breaking = — — (det <1> + det 1 ) + — x tr [<M> 1 ] (det <f> + det 1 ) + — y (det d> + det <I>1) 2 

+w (tr [d^'t 2 d^*t- 2 \ +h.c .). (2.3) 

where the building block <t> = \/2((>o — 4o)4 + V^-iXi + /</>/)£' with to = I /2 and f = a,/2 trans¬ 
forms as 

d> -> e 2i ° A V3>R (L G SU L ( 2), R G SU R ( 2), 0 A G Re), (2.4) 

under chiral and £4(1) transformation. Under the working hypothesis that this system undergoes 
second order phase transition, the existence of an infrared-stable fixed point (IRFP) is examined in 
d = 4 — e dimension with e = 1. 

Rewriting eq. ( |2. 1[ ) in terms of components, (j) a and we obtain 

•^total =(l+w)^(<9 M (/» a ) 2 + ^ (m 2 -y) 0 O 2 + (1-W)^(a^) 2 + ^ (m 2 + y)^ a 2 

+ y [A(^ 2 ) 2 + (A-2x)(^ 2 ) 2 + 2(A+g 2 -z)0« 2 Z6 2 -2g2(te) 2 ] , (2-5) 

1 The early version is Ref. [|6j] 
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where A = gi + £ 2/2 +x + v, z = x + 2y and a runs 0 to 3. In the following, w = 0 is taken for 
simplicity. Non-zero ca breaks degeneracy between <j> a and % a as 

2 2 C A 2 2 , C A ^ 

my, =m ——, m x =m +—. ( 2 . 6 ) 

We take c A > 0 and mp = 0 to realize QCD at T c . Then, Xa has a mass of m x 2 = ca- 

It is naively expected that the massive fields Xa decouple from the system in the IR limit. If it 
is the case, the system will eventually end up with well known 0(4 ) LSM and hence undergo the 
second order phase transition inhering in the 0(4) universality class. 

In order to trace the effect of the massive fields, we take a mass dependent renormalization 
scheme 2 . /3 functions thus obtained are 

j 8 A = R — = - eA + 2A 2 + -/(A)(4A 2 + 6A| 2 + 3 g\ - 8A£- 6g 2 z + 4?), (2.7) 

d On 1^1 ~ 1 /v 

Pg 2 =- £ S 2 + ^^g 2 + -^f(fL)g2(^-2x) + -h(fL)g 2 (4?i+g2-4z), (2.8) 

& = J u^ = -£x + 4f(fi)(Xx-x 2 ) 

+ ^(1 - f(A))(8X 2 - 6A g 2 - 3 g\ + 8Az + 6 g 2 z - 4z 2 ), (2.9) 

fa = = - £z+ ^(2A 2 - Xg 2 + 2A£) - ^/t(A)(4A 2 - 3g 2 - 8Az + 4£ 2 ) 

djl 2 6 

+ \f( A) (-2A 2 + 3Ag2 + 3 g\ - 2A£ - 6 g 2 z + l2Xx + 6g 2 * - 12x£ + 4z 2 ) (2.10) 
6 

where the couplings witlTdenote dimensionless couplings, /./ = Li /cj is the renormalization scale, 
and 

/(A) = 1 -^7TTF" c,an 73TF ACA> = 1 “^ los[I+A21 ’ <211) 

with /(A —» 0) = A 2 /3, h (A -> 0) = A 2 / 2 , and /(A —»°°) = h(p. -> °°) = 1. 

The resulting RG flows can be classified into two cases, depending on the asymptotic behavior 
of A in the —$■ 0 limit. One case tends to be realized when ca is small, where all the couplings 
diverge. Thus, in this case the transition will be of first order. The other case tends to happen 
when ca is relatively large. In this case, all the couplings except for A diverge, but A approaches 
Airfp = e/2 that coincides the IRFP of the 0(4) LSM. 

The asymptotic behaviors of the diverging couplings are calculable and turns out to be 

g2,asym(jU) = Hm | 2 (m) = (2.12) 

,u-t0 

X a sym (jl ) = = ^glasym (M), (2-13) 

2 For two-point functions, the on-shell scheme is applied, and hence c A does not run. 
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Zasym(/-0 — ^§2,asym (M ) > 


(2.14) 


where c is a dimensionless constant depending on the initial conditions of the RG flow. 

With above results, we can calculate the approaching rate of X to its IRFR Towards the infrared 
limit (/t —7 0), becomes 


lim n % =-eA+2A +/(A) ^§2,asym 82 , 


1 


ji->0 


= — eX + 2% 2 + 

24 

Then, the approaching rate is found to be 


^asym^asym 




asym 


+ ... 


lim X(n) — Airfp 


1/3 


Usually, the approaching rate is discussed through the derivative of the /3 function, 

d h\ 


CO = 


dX 


1 A=Airfp’ 


(2.15) 


(2.16) 


(2.17) 


which is one of the universal exponents 3 , and the 0(4) LSM and the (4(1) broken LSM take 


®0( 4) ^ a ^4 ^{/^(ljbroken 1 Xe/3, 


(2.18) 


respectively. 

General argument of renormalization group tells that co is determined by the RG dimension 
of the leading irrelevant operator. In 0(4) LSM, (0 n 2 ) 2 is the one. In the (4(1) broken LSM, 
the leading irrelevant operator is not evident but should not be the same as the 0(4) case because 
(0 0 ( 4 ) / 0 ) !Ja (i )broken ■ To identify the leading irrelevant operator, we calculate co with §2 = 0 as 
a trial, which corresponds to omitting (fl>aXa) 2 [see eq.(2.5)]. Then, co = £ is obtained. From this 


result, it is concluded that the operator (faXa) 2 is the leading irrelevant operator in the (4(1) broken 
LSM. Notice that this operator is not invariant under the 0(4) rotation in the 0(4) LSM, which acts 
only on c/) a but leaves Xa unchanged. 

The above discussion may sound strange because, if the decoupling theorem [|8, j9j] holds in 
this system, the massive field Xa should not affect the low energy behavior of the system. However, 
the decoupling theorem becomes non-trivial when a theory contains a dimensionful coupling 101- 

Si nee one of the universal exponents, CO, turns out to be different from that in the standard 
0(4) LSM, it is clearly interesting to see the other exponents (v and t] in the standard notation). 
We have calculated these exponents in the (4(1) broken LSM through one-loop, and found that 
they are the same as those in the 0(4) LSM. Two-loop calculation is ongoing. 

Whether X flows into 4 rfp or not depends on not only c,\ but also the initial conditions for the 
couplings. We have determined for various values of ca the parameter space of the initial condition 
(attractive basin) with which X —> Airfp. It is found that the attractive basin shrinks as ca decreases. 

The scheme dependence of the low energy behavior was also examined by taking the MS 
scheme as an alternative. Since the MS scheme does not take care of the finite mass of Xa a nd 


3 For exceptions, see Ref. [7]. 
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hence the /3 function does not feel the decoupling, we compared the low energy behavior of the 
four-point Green’s function of (j) a between different schemes. No scheme dependence is observed 
at the one-loop level. For details of the above considerations, see Ref. [jS]. 


3. Effective potential 


In order to gain more understanding, we are calculating the effective potential of the C/ A (l) 
broken LSM. Our setup indicated by the working hypothesis and m<j, = 0 corresponds to the system 
slightly above T c . Thus the order parameter (<f>) must be zero. By looking at the effective potential, 
we check whether this is indeed the case. In the following, we concentrate on the case where 
(Xa) = 0 for simplicity. The calculation of the general case is not completed yet. We choose 
the vacuum expectation value as ((j) a ) = {0 C /, 0,0,0}. The explicit one-loop calculation yields the 
effective potential, 

Vcs(M/n e = yA(At)^ + V^ loo ^</> c /)/At e + ^- loop (0 c/ )/At^ (3.1) 


where the first term represents the tree level result, and the second and third terms are the one- 
loop contributions from the diagrams containing (j) a and Xa i n the loop, respectively. Neglecting 
Vy lo ° p (0 c/ ), V eff agrees with that of 0(4) LSM. The explicit expressions for the one-loop contri¬ 
butions are 


v'T loop (fa)/M' 

'T loop (fa)/M 8 


n- 


A 2 (a 0 log 


4 


tycl 1 




+ I + £.lo g 3 ) (j) cl \ 


(3.2) 




-2e 


4(47t) 2 


(«h) log[mj/c A ] + 3 (m 2 ) log[m 2 /c A ] 

- ( ~z(fl)) 2 + 3(X(fl) + gi(fl) ~z(fl)) 2 } 


dt, log 




CA 


+ \) <t>d 4 




C3.3) 


where 

m? = ^ 2 /i e (A(/i)-£(/i))0 d 2 + c A , ml = ^n 2 fl e (X(fi)+g2(fl)-z(fl))^ c 2 + c A . (3.4) 
By solving the RG equation 





dgi 



V'eff(0c/) = 0, 


(3.5) 


we implement the RG improvement, which is realized by simply setting /./ = <p (; / 2 ^ 2 in the 
unimproved expressions. Where {g,} = {X,g2,x,z}- Because there are two scales <j> c i and c A , we 
need more powerful technique for the resunrmation. Owing to the RG improvement, we can easily 


5 














Renormalization group flow of linear sigma model with UA( 1) anomaly 


Tomomi Sato 


c a /A 2 =1 ,(A,(A),g 2 (A),x(A),z(A))=(0.05,0.5,0,0) 



Figure 1: Effective potential and contributions from each term. Parameters are set to be ca/ A 2 = 1 and 
(A (A), § 2 (A), x(A),|(A)) = (0.05,0.5,0,0). This is inside of the attractive basin, e is set to 1. There is a 
non-trivial minimum produced by higher order contribution in v\ loop . 


see that the effect of V)' loop is just a shift of the coefficient of the tree potential. In the attractive 
basin, A (ft) goes to the positive constant £/2 in jU —^ 0. Therefore, vj loop does not produce a 
non-trivial m in imum as seen in the numerical calculation. 

Fig. [l] shows shapes of the RG improved effective potential drawn by numerical calculation. 
The largest contribution near the origin comes from V% loop , and it produces the non-trivial min¬ 
imum. In this region, the loop contribution is larger than the tree contribution, the perturbative 
calculation is disputable. To get more information about the minimum, we estimate the behavior 
of V% loop near the origin. Since gz^r/ 2 and z^d 2 converge to zero in /./ = </) r / 2,/ ' 2 ~ e) —y 0, we can 
expand V% loop around (j) r i = 0 as 


^■ l00P (fe)/M e «^M e {(A-z) 3 + 3(A+| 2 -z) 3 }^- 
y Ca 

~ C{(A -z) 2 + 3(A + 82 -z) 2 }fe 4 - + • 
6 j c A 


(3.6) 


where the first term comes from (m 2 ) 2 log [nv\/cA and (m 2 ) 2 log \tn\lcp\ in eq.(T3). These terms 
have lineal - and quadratic in couplings , though they are canceled by other terms. The second 
term comes from the one including log[l + q (I — q )/j 2 /ca]. Using the IR asymptotic behaviors of 
couplings (eq.([2.12|)-(|2.14])), 


limu e {(A-z) 3 + 3(A+| 2 —z) 3 }—=-^c 3 c] e V = -^c 3 cf ‘<ju 6 3e , 




ca 

2i j, 4 ^ 5 2 „|e-l ,,2 -Iea 4 


lim{(A - z) 2 + 3(A +| 2 - =^c z c A 

p->o ca o 


0d 4 = |c 2 c* e V c / 8 
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In the rightmost in each line, we use fl = <ft r 2 ^ 2 e '. The first term of cq([3. 6j) is proportional to (j ) r 3 
and the second term is proportional to (j) c i +1 in £ —>1, The first term gives larger contribution than 
( p c / 4 . This term becomes negative when c > 0, therefore another minimum other than (j) r / = 0 exists 
as shown in the numerical calculation. However, this contribution is cubic in couplings. Because 
we carry out the calculations to quadratic order, the cubic order is subject to unknown two-loop 
coefficients. To know whether the non-trivial minimum arises or not, the two-loops calculation is 
needed. 

4. Summary 

In the £ expansion, we calculated the RG flow and the effective potential of the £4( 1) broken 
LSM. There are two cases of the flow depending on size of the £4(1) breaking. With small £4(1) 
breaking, the flow tends to diverge in the IR limit. Increasing the breaking, one of the coupling 
converges to the same value with the IR fixed point of the 0(4) LSM. In this case, the critical 
exponents, co, V, t] are calculated to one-loop. 0) turns out to differ from that in 0(4), while 
the others are found to be the same to one-loop order. This difference comes from the leading 
irrelevant operator in the £4(1) broken LSM and can be interpreted as the non-decoupling effect. 
We determined the attractive basin in parameter space of the £4(1) broken LSM. This attractive 
basin shrinks as the size of the £4(1) breaking decreasing. There is a non-trivial minimum in the 
RG improved effective potential. However, this minimum is produced by the order in couplings 
higher than we have performed. We need to cany out the higher loops calculation to establish the 
location of the minimum. 
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